Viscosity and Thermal Relaxation for a resonantly interacting Fermi gas 
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The viscous and thermal relaxation rates of an interacting fermion gas are calculated as functions 
of temperature and scattering length, using a many-body scattering matrix which incorporates 
medium effects due to Fermi blocking of intermediate states. These effects are demonstrated to 
be large close to the transition temperature T c to the superfluid state. For a homogeneous gas in 
the unitarity limit, the relaxation rates are increased by nearly an order of magnitude compared 
to their value obtained in the absence of medium effects due to the Cooper instability at T c . For 
trapped gases the corresponding ratio is found to be about three due to the averaging over the 
inhomogeneous density distribution. The effect of superfluidity below T c is considered to leading 
order in the ratio between the energy gap and the transition temperature. 
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I. INTRODUCTION 

The properties of resonantly interacting fermions have 
been under intense investigation during the last few 
years, stimulated by the ability to manipulate the in- 
teratomic interaction through the use of Feshbach res- 
onances. In this way it has been possible to create 
molecular Bose-Einstein condensates from weakly bound 
molecules consisting of two fermionic atoms. The molec- 
ular condensates are formed when the scattering length is 
positive, corresponding to an effective repulsion between 
particles. If the interaction is attractive (the scatter- 
ing length negative) the atoms may pair up in a man- 
ner similar to the way in which electrons form Cooper 
pairs in superconducting metals. The important dimen- 
sionless parameter for an interacting, degenerate fermion 
gas is the product of the Fermi wave number fcp , given 
by fcp = 37r 2 n where n is the fermion density, and the 
scattering length a. For attractive interactions, and for 
k-p\a\ much less than unity (which we refer to as the 
weak-coupling limit), the transition temperature T c to 
the superfluid phase is very much less than the Fermi 
temperature Tp . Near Feshbach resonances it is possible 
to achieve values of fcp|a| that are much larger than unity 
(the strong-coupling limit), and in this case the expected 
transition temperature is comparable to the Fermi tem- 
perature, typically about one quarter of Tp, depending 
on the model used for the effective scattering matrix. 

In this work we demonstrate how many-body effects 
have a significant influence on the viscous and thermal 
relaxation rates at temperatures near the transition tem- 
perature to the superfluid state. The many-body effects 
are due to Fermi blocking of intermediate states in the 
effective scattering matrix. Such effects arc known to in- 
fluence equilibrium quantities such as chemical potential 
and total energy, but - as we shall see below - their in- 
fluence on properties like thermal or viscous relaxation is 



much more pronounced. Thermal relaxation rates may 
be determined directly experimentally Q, Q , and viscous 
relaxation rates can be extracted from the attenuation 
of collective modes. Our work is an extension of Q, in 
which the viscous relaxation rate was calculated using 
the vacuum scattering matrix, and we shall refer to this 
paper as I in the following. 



II. VISCOSITY AND THERMAL RELAXATION 

In this section we discuss the relation between the ther- 
mal and viscous relaxation rates for a two-component 
Fermi gas trapped in a harmonic-oscillator potential, 
with an equal number of particles in each component. 
The interaction between the fermions is taken to be at- 
tractive, a < 0, and we consider fermions in their nor- 
mal state, at T > T c . For bosons in a harmonic trap, at 
temperatures above the Bose-Einstein transition temper- 
ature, the thermal and viscous relaxation rates have been 
shown to be simply related, differing only by an overall 
factor of two |^. As one might expect, the same relation 
applies to the fermion gas, since the factor of two origi- 
nates in the equality of kinetic and potential energies for 
a particle in a harmonic trap. We briefly indicate the 
line of reasoning which leads to this result and refer the 
reader to 3 for details. 

We consider a gas of particles with mass m, trapped 
in an axially-symmetric harmonic-oscillator potential of 
the form 



y(r) = ^iri+^z 2 ), 



(1) 



where 



y 2 . For such a system we define the 



spatially averaged viscous relaxation rate by 

1 _ Jd 3 r(X,H[X]) 
r rl jd^r(X,X) ' 



(2) 



where the symbols have the same meaning as in I. In 
particular H is an integral operator derived from the col- 
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lision integral of the Boltzmann equation, and X is a trial 
function given by 



X = VyP X , 



(3) 



with p being the particle momentum and v the group 
velocity, equal to p/m in the normal state. In carrying 
out the spatial average we treat the trapping potential 
in the Thomas-Fermi approximation via a locally vary- 
ing chemical potential /i(r) = fj, — V(r). The time t v 
determines the nature and the damping of the collective 
modes in the trapped atomic gases 

The measurement of thermal relaxation rates has 
proven a useful tool for investigating the properties of 
trapped gases [J [2]. In such an experiment, the gas is 
"heated" preferentially in one spatial direction, and the 
relaxation toward equilibrium with a uniform tempera- 
ture in all directions is determined from the time evolu- 
tion of the rms cloud radii. The thermal relaxation rate 
can be extracted from the time-dependent aspect ratio. 
In analogy with the viscous relaxation rate, one may ob- 
tain a variational expression for the relaxation time of 
temperature anisotropies in a trap Q , 



1 _ J d 3 r{<P T ,H{<S> T }) 



f d 3 r($ T ,(f> T ) 



where the trial function is 



$ T = Pt - P 2 /'S + -^(^iz 2 - wjri). 



(4) 



(5) 



The momentum-dependent part of this trial function in- 
volves only a I = 2 spherical harmonic in momentum 
space. If the spatial part of the trial function were ne- 
glected, the viscous and thermal relaxation rates would 
be identical, since the trial function @ also involves only 
I = 2 spherical harmonics. Since collisions relax only 
momentum anisotropies, the spatial part of © does not 
contribute to the numerator in whereas it gives the 
same contribution to the denominator as the momentum 
part. It therefore follows from comparing @ and (@J 
that the thermal and viscous relaxation rates are related 
by 



T T = 2t„ 



(0) 



The viscous relaxation rate is thus directly accessible ex- 
perimentally via measurements of the thermal relaxation 
rate as in Q, 13 ■ 

III. COLLISIONS 

In order to calculate the viscous relaxation rate, we 
need the particle-particle scattering cross section. The 
differential cross section da/dQ for s wave scattering is 
related to the scattering matrix T by 



da_ m 2 \T\ 2 

dh 



(Anh 



2 \2- 



(7) 



An often used approximation close to a resonance is to 
use the so-called unitarized vacuum scattering matrix 



4irh 2 a 



1 



1 



iqa 



(8) 



where a is the scattering length and Tiq the relative mo- 
mentum of the scattering particles. In the present paper, 
we shall extend the calculation carried out in I of the vis- 
cous relaxation rate to cover the temperature region near 
T c where medium effects absent in the approximation © 
turn out to be important. 

Close to a Feshbach resonance, the scattering length 
can be written in the phenomenological form 



1 



AB 
B-Bn 



(9) 



Here, AB is the width of the resonance, db g the back- 
ground scattering length coming from non-resonant scat- 
tering processes, and Bq the position of the resonance. 
For atomic gases close to a Feshbach resonance, a multi- 
channel effective theory for atom-atom scattering was 
recently presented The scattering matrix including 
medium effects was shown to be given by 



T(u,K) 



bg 



1 



A/jA-B 



hu>+h(u>,K.)-AiJ,(B-B ) 



-r bg n(w,K) 

where KK. is the center-of-mass momentum of the scatter- 
ing particles and uj the frequency. Here Co = lo — UK 2 /4m 
and 7b g = Airabgfi 2 /m, while Afi is the magnetic moment 
of the Feshbach molecule. The pair propagator n is 



n( w ,K) = 



l-/°(±K-q')-/°(±K + q') 



1 



iS 



h 2 K 2 _ H 2 q' 2 



h 2 q' 2 jra — i5 



(11) 



where / is the equilibrium Fermi function. The quan- 
tity h(u>,T£) in l|l(J|) describes effects coming from the 
fact that the Feshbach state is a composite two-fermion 
object, but it can be neglected for most resonances of 
experimental interest and we therefore set it equal to 
zero, h(ui, K) = 0. For the calculation of transport coef- 
ficients such as the viscosity, the scattering matrix has to 
be evaluated on-shell with Huj equal to the energy of the 
two scattering particles, i.e. huj — h 2 K 2 /4m+h 2 q 2 /m. In 
a vacuum, the on-shell version of (|10|l with h(u),K) = 
reduces to 



T 



A/jA-B 



h 2 q 2 /m-An(B-B ) 



(12) 



tab g q 



We shall use (|1(J|) for calculating the low-temperature vis- 
cosity for a resonantly interacting gas. 
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The approximation (JSJ) for the scattering matrix ne- 
glects two effects as compared to the many-body expres- 
sion (|10|l . Firstly, even the vacuum limit (|12Jl of jllOJl is 
not identical to ©; they agree only when the ft, q 2 /m 
term in (|12|l can be neglected. This term corresponds 
to an effective range of the atom-atom interaction. Here 
we consider only broad resonances such as the one at 
B ~ 201. 6G for 40 K or the one at B ~ 830G for 6 Li, for 
which the effective range can be safely neglected |a, Q • 

Secondly, the approximation given by JSJ neglects 
medium effects, which are included in i|10[l via the prop- 
agator ljll|l . These effects come from the occupation of 
open-channel states. One effect of the medium is to shift 
the resonance position away from its vacuum value Bq p| . 
More importantly for the present purpose, medium ef- 
fects cause a significant increase in the scattering rate 
for T — + due to the Fermi blocking of the pair states 
into which the molecular state can decay. The decay of 
the molecular state in a vacuum is the origin of the iqa 
terms in the denominator of (jHJ and (|12|) . These terms 
are replaced by T^ g H in l|lUfl when medium effects are 
taken into account. At low temperatures, this term is 
significantly reduced from the vacuum value due to the 
Fermi-blocking factors in (|llf> : the lifetime of the res- 
onant state is increased leading to a corresponding in- 
crease in the scattering rate. We shall demonstrate that 
this effect is significant close to the superfluid transition 
temperature T c . 



A. Low and high temperature limits 

Before we present any results, let us first discuss what 
can be calculated analytically. For low temperatures 
T/T F < 1 (fcT F = h 2 k 2 ,/2m), the Fermi factors in the 
collision integral (see Eq. (13) in I) restrict the avail- 
able phase for the scattering particles such that their 
momenta are close to the Fermi surface. This simpli- 
fies the integrations and assuming the simple unitarized 
vacuum form J3J) for the scattering matrix (T = Tuni), an 
analytical result for the viscosity can be obtained using 
standard Fermi liquid tricks |g. A lengthy but straight- 
forward analysis yields for fcp|et| <§; 1 and T — > 



3nh 



V = 



8n(k F a) 2 (kTf 



(13) 



with e F = h k F /2m. Equation @3J| is based on T < T F 
and the assumption T = T un i- This assumption works 
well for weak coupling with T c <C Tp: there is then a 
temperature regime for which T <C Tp such that all scat- 
tering momenta are confined to the Fermi surface, and 
T ^ T c such that medium effects can be neglected. How- 
ever, for strong coupling T c becomes comparable to T-p 
and there is no temperature regime for which T -C Tp 
and T > T c . 

In the classical limit T 3> Tp, an analytical expression 
for the viscosity can also be obtained using T = T un i\ 



5\/ mkT 
320? 



5 V 'irmkT 

a- 2 , T«T a 

3mkT/h 2 , T>T a 



(14) 



IV. HOMOGENEOUS SYSTEM 



with kT a = Ti 2 / (ma 2 ) and 



Ana 2 



We now present results for the viscosity of a homoge- 
neous system. The viscosity is calculated as a function of 
temperature for a constant density n. We calculate the 
thermodynamic potential f2 for fixed n within the lad- 
der approximation as a function of temperature, taking 
into account two-body scattering processes with a scat- 
tering matrix given by l|l(J|) . The chemical potential fi 
is determined by the condition —d^U, — n. Our calcu- 
lation of £1 is identical to the one described in Ref. 
where more details are given. Once the chemical poten- 
tial is determined, we calculate 77 using Eq. (16) of I with 
the cross section given by Q and (|l(Jfl . The calcula- 
tion of the chemical potential and the viscosity requires 
multi-dimensional integrations which are performed nu- 
merically using two different methods to check the ac- 
curacy. The first method uses convoluted ID integra- 
tion routines and the other is based on multidimensional 
Monte-Carlo integration. Both yield the same result for 
the chemical potential and the viscosity within less than 
1%, thus confirming the accuracy of the numerical calcu- 
lations. 



dxx'e- x (l+x 2 T/T a y 



(15) 



the momentum averaged cross section. As opposed to 
the low-temperature limit, the assumption T = T unl is 
reasonable for T Tp, even in the strong-coupling case, 
since medium effect on the scattering matrix are sup- 
pressed at high temperatures. 



B. Numerical results 

Fig. H shows the viscosity r\ in the weak-coupling limit 
fcp|a| <C 1. We have chosen parameters in (|1U|) corre- 
sponding to a generic weak interaction characterized by 
background scattering with k F ah g — —0.1, AB/(B — 
Bo) = 0.0012, and a negligibly small effective range. 
For weak interaction, the scattering matrix reduces to 
Airh 2 ah s /rn and there are no effects of the medium as 
long as T is well above T c <§C Tp. In fact, the solid line 
in Fig. n is indistinguishable from the result obtained if 
we had used 4irh 2 ahg/m for the scattering matrix of i|10|) 
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ical BEC-BCS crossover results based on a similar ap- 
proximation for the thermodynamic potential [{J, Il0| . 

Just above the critical temperature T c , the viscosity is 
reduced almost by an order of magnitude (a factor w 7.5) 
as compared to the value obtained using the vacuum scat- 
tering matrix. The reduction reflects the increase in the 
scattering rate due to Fermi blocking of the molecule de- 
cay as explained in Sec. 1 1 1 1 1 These medium effects cause 
the zero center-of-mass (COM) momentum T matrix to 
diverge at T c , which is the signature of the Cooper in- 
stability. Thus, we expect the increase in the scattering 
rate due to the medium to be significant for tempera- 
tures close to T c , regardless of the value of the coupling- 
strength. This is confirmed by numerical calculations for 
different values of the coupling: Medium effects on the 
scattering rate and thus the viscosity are significant when 



FIG. 1: (Color online) The viscosity r\ in units of nh for back- 
ground scattering. The solid line is the numerical result with 
the scattering matrix ill (Ifi . the dotted lines are the low and 
high temperature limits. 
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T/T 

c 

FIG. 2: (Color online) The viscosity r\ in units of nh for reso- 
nant scattering as a function of temperature for T > T c . The 
solid line is the numerical result with the full scattering matrix 
(1 1 the dashed line the numerical result using the vacuum 
approximation ©. The dotted line is the low-temperature 
result (1141 . The inset shows the chemical potential. 



as expected. We see that the viscosity approaches the 
low- and high-temperature forms given by (|13|l and i|14|) 
as expected, with a minimum located at T w 0.7Tp. 

Fig. [21 shows the viscosity i] in the strong-coupling 
limit. We have chosen parameters in (|10|) correspond- 
ing to a generic resonant interaction characterized by 
hp a = —If. 8 and a negligible effective range. For this 
set of parameters very close to resonance, we find the 
critical temperature for the superfluid transition to be 
T c w 0.26Tf in excellent agreement with earlier numer- 



In the high temperature limit on the other hand, the 
medium effects on the scattering are suppressed as ex- 
pected and the viscosity approaches its classical value as 
can be seen from Fig. [21 We conclude that the medium 
increases the scattering rate significantly as compared to 
the vacuum value for temperatures close to T c due to the 
Cooper instability. 

Note that since kT c / n(T c ) w 0.6 is rather large for 
strong coupling, the T -2 increase in the viscosity coming 
from the phase space blocking in the collision integral is 
not observed for T > T c in Fig. 

The above results show that the minimum viscosity 
for a Fermi gas in the normal phase for k-p\a\ — > oo and 
T = T c is 

?7min = anh (17) 

with a ss 0.15. This minimum is significantly less than 
the value a ~ 1.1 obtained when the vacuum scatter- 
ing matrix given by © is used. We caution that the 
minimum value of rj/nh obtained in (|17|l depends on the 
validity of our starting point, the semiclassical kinetic 
equation. A simple relaxation-time approximation to the 
collision integral, yielding the low-temperature viscosity 
7] mnvpT/5, would result in i] w OAnh if r = %/e-p- 
However, for such small values of r the semiclassical ki- 
netic equation may need modification due to the spectral 
broadening of the momentum states. It is interesting to 
note that a related problem, the viscosity of a hot quark- 
gluon plasma, was recently analyzed within the Kubo 
formalism [flj which yielded results that differed signifi- 
cantly from those obtained from the Boltzmann equation. 

C. Effects of superfluidity 

The viscosity of a superfluid gas of fermions has 
been investigated extensively in the context of superfluid 
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3 He Below the transition temperature T c , the vis- 

cosity drops rapidly as a function of temperature, the 
relative decrease in viscosity being of order A/kT c . The 
physical reason for the large drop is the change in the 
quasiparticle dispersion relation and the modification of 
the collision operator. For the purpose of illustration 
we shall here discuss the two effects separately, although 
they must of course be treated together in a consistent 
calculation. 

Firstly, the group velocity is changed in the superfluid 
state, since the dispersion relation of Bogoliubov quasi- 
particles is 



e = + a 2 , 



(18) 



where £ = p 2 /2m — /i and A is the energy gap. The group 
velocity v s = dE/dp in the superfluid state is thus 



V = p«F 



i 

E' 



(19) 



where p is a unit vector along p and vf = Tik-F /m is the 
Fermi velocity. This affects the velocity component v y 
appearing on the left hand side of the kinetic equation 
and in the expression for the momentum current density. 
To illustrate the effect of the change in group velocity, let 
us make a relaxation time approximation to the collision 
operator occurring in the Boltzmann equation. The shear 
viscosity is obtained by assuming a flow velocity u of 
the form u = (u x (y), 0, 0) and linearizing in the gradient 
diix/dy, as a result of which the kinetic equation for the 
non-equilibrium distribution function / becomes 



du x 



dy VyPx de 



df° f-f 



(20) 



The shear viscosity rj relates the momentum current den- 
sity U xy , given by 



dp 



■VyPxf, 



(21) 



to the gradient of the flow velocity according to H xy = 
—r)du x /dy. Inserting the resulting distribution function 
in (I21|) and carrying out the integration over angles we 
can perform the energy integration analytically, using the 
fact that 



£ 2 



1 - ttA<J(0. 



(22) 



Here the representation of the Lorentzian by a delta- 
function is valid as long as the remaining part of the 
integrand varies slowly on the scale of A, that is pro- 
vided A < kT. Since the value of /°(1 - /°) at the 
chemical potential is 1/4 for a degenerate fermion gas we 
obtain 



V = Vn 



7T A 

1 ~ 1W r . 



(23) 



to lowest order in A/kT c . The quantity rj n = mnvpT/5 
is the viscosity in the normal state at T c . 

Secondly, the probability of two-quasiparticle scatter- 
ing processes is modified by the coherence factors appear- 
ing in the Bogoliubov transformation. Since the number 
of quasiparticles is not a conserved quantity in the su- 
perfluid state, it is also necessary to include processes in 
which one quasiparticle decays into three as well as the 
inverse processes. 

Taking both these effects into account, it is possible 
near the transition temperature to express the relative 
decrease in viscosity in terms of normal-state quantities 
and relate the coefficient in front of A/kT c to the ratio 
between the viscous relaxation time t v and the lifetime 
t(0) of a quasiparticle in the normal state at the Fermi 
energy. One finds [13j that the change in viscosity is 
given by 



r, 4 V 12 t(0) 



t(0) A 
TV) kT r . 



(24) 



The coefficient of the A / kT c term may thus be expressed 
solely in terms of normal-state quantities. For the trial 
function © used here, the ratio between the lifetime and 
the viscous relaxation time is [HI 



r(0) 



x 5 da 
dx — — — 

o VT^x 1 dn 



1 J x d * 1 

dx — — — 

(25) 

The medium cross section diverges at zero COM momen- 
tum at T c . This leads to a non-integrable pole at x = 1 
in both integrals in l|25ll and therefore 



r(0) 



(26) 



This is a universal result, independent of the coupling 
strength, coming from the Cooper instability at T c . 

The viscous relaxation time is obtained from the vis- 
cosity by removing the factor (X, X), where X is the 
trial function proportional to v y p x . By the same reason- 
ing that led to 12: ill it follows that the value (X,X) S of 
the scalar product in the superfluid state is related to its 
value (X,X) n in the normal state by 



(X,X) s ~(X,X) n 1 



7T A 



(27) 



to leading order in A/kT c . Consequently, the change in 
r r) is given by 



St v A 
—1 = 0.0643 — 



(28) 



Due to the smallness of the numerical factor in l|28|l we 
conclude that there is very little change in the viscous 
relaxation rate upon entering the superfluid phase. 
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V. TRAPPED GAS 

We now present results for the thermal relaxation rate 
as a function of temperature both in the weak and the res- 
onant coupling limit for a trapped two-component Fermi 
gas. To keep the number of trapped particles N = 2N a 
constant (there are N a atoms in each hyperfine state), 
we calculate fl in the ladder approx. as in Sec. and 
the chemical potential is determined by the condition 
N = —d^Cl. The relaxation rate is then calculated 
from Q. The presence of the trapping potential V(r) 
is treated in the Thomas-Fermi approximation by a local 
chemical potential fi(r) — fi — V(r). This requires an 
additional integration over the spatial coordinates r, and 
the resulting multidimensional integrals are evaluated us- 
ing a Monte-Carlo routine. 




o 1 1 1 1 1 — 

1 2 3 4 5 6 7 

T/T 



A. Low and high temperature limits 

As for the homogeneous case, the integrals in Q for 
the thermal relaxation rate can be calculated analytically 
in the low and high temperature limits using the unita- 
rized vacuum form T un i for the T matrix. For T/Tp <C 1 
and weak coupling /cf |«| <C 1, we obtain 



1 



16n(kT) 2 , 

= — ma . 

15h 3 



(29) 



For high temperatures, T 3> Tp, we obtain 
1 N muj 3 a 



t t IOtt 2 kT 



2NCu- 



~,3 



2 , T<T a 

5irkT ~ \ h 2 /(3mkT) , T > T a 



(30) 



with uj = uj x uj y uo z and a given by (|15|l . As for the homo- 
geneous case, T = 7^ n j is a good approximation at low 
temperatures for weak coupling only, whereas it is valid 
at high temperatures for all coupling strengths. 



B. Numerical results 

For weak coupling kp\a\ -C 1, the numerical calcula- 
tions confirm as expected that there are no effects of the 
medium (except for very small T close to T c ) and the re- 
sults for the thermal relaxation rate are identical to the 
ones reported in Fig. 2 of I. 

We therefore concentrate on the strongly interacting 
limit. Fig. [3| shows the thermal relaxation rate 1 
as a function of temperature for a strongly interacting 
system with Arp|a| 3> 1. The rate is in units of ep/h 
where ep = hd!(3N) 1 ^ 3 is the Fermi energy for N non- 
interacting trapped particles. We have chosen parame- 
ters in (|10|) corresponding to a resonant interaction with 
kpa = —11.8 and a small effective range fcF^cff < 1. As 
for the homogeneous case, the strong interaction gives 



FIG. 3: (Color online) The thermal relaxation rate in a trap 
in units of ep /% for the resonant coupling case fcp|a| 3> 1. 
The solid line is the numerical result with the full scattering 
matrix in HUH , the dashed line the numerical result using the 
vacuum approximation ||SJ for the scattering matrix. The 
curves are given for T > T c ~ 0.3Tf. 



rise to a superfluid transition below a critical tempera- 
ture T c . We find T c ps 0.3Tf at resonance faki 3> 1) in 
excellent agreement with the results in Ref. [l5| . 

We see from Fig. [3] by comparing the solid and the 
dashed lines, that the medium increases the relaxation 
rate significantly for low temperatures as compared to 
the vacuum prediction: close to T c the rate is approx- 
imately 3.6 times higher than the vacuum result. The 
increase is smaller than for the homogeneous case since 
the effects of the medium are reduced near the edge of the 
cloud where the density is reduced. The reason for the 
significant increase of the rate is of course identical to the 
effect discussed for the homogeneous case: the medium 
increase of the resonant scattering rate for low tempera- 
tures (T — T c )/T c < 1 signals the Cooper instability at 
T c . For high temperatures, the medium effects are negli- 
gible. As for the homogeneous case, the T 2 dependence 
of the rate coming from the Fermi blocking of the avail- 
able scattering states for T/Tp <C 1 is not observed for 
T >T C since kT c /fi(T c ) « 0.5 is rather high. 

Fig. |3| shows the maximum relaxation rate for a two- 
component trapped Fermi gas for fcp|a| — > oo. The large 
difference between the relaxation rate for a weakly inter- 
acting system and a resonantly interacting system should 
be contrasted with the much smaller effects of the reso- 
nant interaction on thermodynamic properties. 



C. Experiments 

Let us comment on the possible experimental observa- 
tion of the effects described above. The thermal relax- 
ation rate can be measured by "heating" the gas pref- 
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erentially in one spatial direction followed by observing 
the time evolution towards equilibrium with a uniform 
temperature in all directions 0, 0] • These experiments 
have already demonstrated the saturation of the scatter- 
ing rate at resonance in the classical regime. By per- 
forming such experiments for decreasing temperatures, 
one should observe a significant increase in the rether- 
malization rate as compared to the vacuum rate when 
(T - T c )/T c < 1. 

VI. CONCLUSION 

Using a variational approach, we have calculated the 
viscosity and the thermal relaxation rate of an interact- 



ing Fermi gas for weak and for strong coupling. Both 
homogeneous and trapped systems were considered and 
medium effects on the scattering properties were taken 
into account. The rates were shown to be significantly in- 
creased for temperatures close to the superfluid transition 
temperature. This is due to the presence of the medium 
and signals the onset of Cooper pairing at T c . We fur- 
thermore considered the effects of superfluidity below T c 
to leading order in the energy gap. The effects described 
in this paper can be measured directly by rethcrmaliza- 
tion experiments. 
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